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Inga Kamp and Gerd-Jan van Zadelhoff 



o 
o 

5 



> 
O 

o 
m 
in 
o 



6 



Leiden Observatory, PO Box 9513, 2300 RA Leiden, The Netherlands 
Received December 18, 2000; accepted May 4, 2001 

Abstract. In circumstellar disks or shells it is often assumed that gas and dust temperatures are equal where 
the latter is determined by radiative equilibrium. This paper deals with the question whether this assumption 
is applicable for tenous circumstellar di sks ar ound young A stars. In this paper the thin hydrostatic equilibrium 
models described by Kamp & Bertoldi (200C) are combined with a detailed heating/cooling balance for the gas. 
The most important heating and cooling processes are heating through infrared pumping, heating due to the drift 
velocity of dust grains, and fine structure and molecular line cooling. Throughout the whole disk gas and dust are 
not efficiently coupled by collisions and hence their temperatures are quite different. Most of the gas in the disk 
models considered here stays well below 300 K. In the temperature range below 300 K the gas chemistry is not 
much affected by T sas and therefore the simplifying approximation T sas = Td ust can be used for calculating the 
chemical structure of the disk. Nevertheless the gas temperature is important for the quantitative interpretation 
of observations, like fine structure and molecular lines. 
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1. Introduction 

In the literature the expression "Vega-type" stars is widely 
used to characterize main-sequence stars of any spectral 
type which show an infrared excess. In order to deal 
with a more homogeneous subgroup this paper concen- 
trates on nearby A dwarfs which are surrounded by dust 
disks including the most prominent and well studied stars 
j3 Pictoris and Vega. 

The literature presents a large variety of models for 
these stars ranging from simple spherical shells with con- 
stant density (Chini et al. 



1990) to detailed disk mod- 



els with power law density distributions and several dust 
components (Pantin et al. 1997 ). These models are used 



In order to study the nature of the circumstellar sur- 
roundings and the ongoing physical processes in more de- 
tail Kamp & Bertoldi ( 2000| ) developed disk models com- 
prising the two components gas and dust. These models 
include a realistic treatment of the UV radiation field and 
the chemistry by the use of a chemical reaction network. 
As a first approximation gas and dust are assumed to be 
effectively coupled by collisions and their temperature is 
derived from the radiation equilibrium of the dust. The 
basic result of that work is an explanation of the non- 
detection of CO radio lines without assuming a general 
gas depletion in the disk. The recent detection of H2 in 
the disk around /3 Pictoris (Thi et al. 2001) can be ex- 



to derive constraints on the properties of the dust compo- 
nents like grain sizes, composition and dust mass in these 
disks. 

Several attempts to detect CO or other molecules 
at radio wavelength in the disks around A stars failed 
(Yamashita et al. |l993j Dent et al. |l995| ; Liseau |1999| ). 



On the other hand circumstellar gas is observed in the 
visible spectra of several A-type dwarfs (e.g. Holweger 



& Rentzsch-Holm 


1995 


; Jolly et al. 


1998 


; Welsh et al 


199S; Holweger et al. 


L999J Roberge et al. 200C) proving 



without doubt the existence of atomic and molecular gas 
(/? Pictoris: CO) in their surrounding. 



plained in the context of the above described chemical 
disk models. 

This paper aims at a more realistic determination of 
the gas temperature in the above cited disk models by the 
use of a detailed heating/cooling balance. The assumption 
of T gas = Tdust will be evaluated. The gas temperature is 
important for the interpretation of observational data via 
modelling of the line emission from the gas phase, e.g. CO 
radio lines, [Cll] and [Oi] fine structure lines, and nar- 
row circumstellar absorption lines on top of the broadened 
stellar profile like the ones seen in the Can K line. 
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2. The circumstellar disk model 



2.1.1. Heating processes 



The disk model and the chemical network is described in 



detail by Kamp & Bertoldi ( 2000 ). We restrict ourselves 
here to a short summary. 

The density distribution in the disk follows from a thin 
hydrostatic equilibrium model 



n(r, z) = rii (r/Ri) 



-2.5 -z 2 /2h 2 



(1) 



Here the dimensionless scaleheight H = h/r is assumed to 
be 0.15. The inner radius of the disk Ri is fixed to 40 AU, 
the outer radius R to 500 AU. The dust temperature fol- 
lows from radiative equilibrium assuming large spherical 
black body grains of size a 



[ dust 



= 282.5 (LJL Q ) 1/5 (r/AU)~ 2/5 (a/fimy 1/5 , (2) 



with the stellar luminosity in units of the solar one L^/Lq. 
Such a disk model leads to a surface density which follows 
a r~ 1,2 power law in reasonable agreement with the "lit- 
erature exponents" of brightness profiles (Hayashi et al. 
1985| ; Dutrey et al. |1996| ; Augereau et al. |2001| ) ranging 
from —1 to —1.5. The radiation field F v is described by 



an ATLAS9 photospheric model (Kurucz 1992 ) for the 
appropriate stellar parameters. 

We derive a stationary solution for the chemistry us- 
ing a chemical network, which consists of 47 atomic, ionic 
and molecular species that are related through 260 gas- 
phase chemical and photoreactions. A number of reactions 
is treated in more detail like H2 and CO photodissociation, 
and C ionisation. The only surface reactions incorporated 
are H2 formation and freezing out of CO on cold dust grain 
surfaces (see Kamp & Bertoldi 200C for further details). 



2.1. The gas temperature 

Kamp & Bertoldi (|200C|) assumed that gas and dust are 
effectively coupled by collisions and hence both have equal 
temperatures, determined by the radiative equilibrium of 
the dust. This assumption will be critically evaluated by 
implementing a detailed heating-cooling balance for the 
gas phase. The particle densities depend on the gas tem- 
perature and vice versa. Hence the energy balance of the 
gas and the chemical network have to be solved simulta- 
neously. The gas temperature is determined by a detailed 
energy balance Y = A, where V and A are the sum of 
all relevant heating and cooling rates respectively. To find 
the solution to T — A = 0, we use Ridder's method, a 



root bracketing algorithm (Press et al. 1997), in a slightly 
modified way: as starting values for the gas temperature, 
we use twice and half the value at the previous radial 
point. Each change of temperature due to the root brack- 
eting is immediately followed by a solution of the chemical 
network. If the root bracketing algorithm gets stuck, the 
routine starts to subsequently extend again the bracketing 
interval. In smooth areas convergence is normally achieved 
after 5 iterations. 



In the following we give a short description of the heating 
processes taken into account for the determination of the 
gas temperature. 

Photoelectric heating. The photoelectric emission from 
grain surfaces is a major heating source for the gas 
(Watson [1972 ; Draine 1978). Dust particles in the cir- 
cumstellar disks considered here have typical sizes of a 
few microns; hence the properties of these large grains — 
large compared to the wavelength of UV radiation — are 
well described by those of bulk mate rial. This approach 
differs from the one of Draine (197S ) in that we deal in 
this paper with micron sized dust instead of small inter- 
stellar dust particles. Two different grain compositions are 
studied: graphite and silicate. 

The work functions of these materials are taken from 



the experimental work of Feuerbacher & Fitton (1972), 
w = 4.7 eV for graphite and w — 8.0 eV for silicate. 
Analytical fits to their data describe the photoelectric 
yields Y as functions of the energy of the incoming UV 
photon hv 



\ogY{hv) = < 



' -19.0 exp(-0.263/«/) - 1.44 

for graphite 

-300.0 exp(-0.577/u/) - 1.06 

for silicate . 



(3) 



The kinetic energy spectrum of the emitted p hotoel ectron 
is crudely approximated by Eq.(3) of Draine ( 1978 ) 



f(E,hu) = (hv-wY 



if < E <hv 



w 



(4) 



Here E denotes the kinetic energy of the escaping photo- 
electron. 

The photoelectric heating rate depends on the energy 
of the impinging UV photon and on the grain potential 
U. The currents that determine the grain potential are 
impinging electrons and UV photons. The contribution of 
impinging ions, mainly protons, is negligible, because hy- 
drogen is mostly neutral or even molecular in these disks. 



The incident electron current is (Draine 1978) 



/ kT 
n e s e l (•■ 

\2irm c 



if 6 > 



if 4> < , 



(5) 



with (j> = eU /kT, m e the electron mass and n c the electron 
density. The sticking probability of electrons to grains s e is 
assum ed to be 0.5 following the argumentation of Draine 
( |1978| ). 

Photoelectrons will be emitted from th e grai ns at a 
rate per grain surface area given by (Draine 1978 ) 



hbsY(hv)f(E, hv)F v dv\dE , (6) 
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Fig. 1. A comparison between the exact photoelectric 
heating rate (plus signs), the analytical approxima- 
tion for large graphite grains (Eq.(||): solid line, a = 
2.34 1 - 23 cm 2 H-atom -1 ) and the ana lytica l approxima- 
tion for small grains (Bakes & Tielens [1994 dashed line) 
as a function of the grain charge parameter assuming ISM 
conditions, x = 1) an d two different gas temperatures 
T = 100 K and 10 4 K. 



with E min = for U < 0, E min = eU for U > and 
-E-max = tomax — tu. The absorption efficiency of the dust 
grains Q a bs is assumed to be 1 independent of wavelength 
in the UV. The frequency vth denotes that of impinging 
UV photons at the threshold w. The equilibrium grain 
potential U can be derived from equating (|5|) and (^|) . 

Given the grain potential and the stellar UV radia- 
tion field, the photoelectric heating rate can be determined 
from 



Ti = 4nHC 



Bmin 



5 abs 



Y{hv)f{E,hv)F v dv 



(E - eU) dE erg cmT A s" 1 , (7) 



where «h and a denote the particle density of hydrogen 
and the grain geometric cross section per H nucleus. 

In order to make the calculation of the photoelectric 
heating rate efficient enough for the inclusion into the 
heating/cooling balance, the rate is approximated using 
a function similar to Bakes & Tielens (1994) 



rv 



10 a e x nn 



(8) 



gas 



< X < 10 , where the flux is measured in units 



of the Habing fll968| ) field between 912-1110 A (F H = 



)■ 

The photoelectric efficiency e for graphite is approxi- 
mated by 



e = 



610" 



1 + 2.510- 4 



,0.95 



y (i(r 4 r) 
l + io- 2 x 



(9) 
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Fig. 2. Same as Fig. 0but for large silicate grains 

where y depends on the so-called grain charge parameter 
x = x T ' 5 /n c 

if x < 10" 4 
if IO" 4 < x < 10- 1 
if IO" 1 < x < 10 2 
if x > 10 2 




(10) 



This formula yields photoelectric heating rates accurate 
to within 10-25%. 

The photoelectric heating rate for silicate grains is 
slightly different from the one for graphite and reads as 
follows 



r po = 2.5 10~ (T e x «h 
with 

610~ 2 

e = 



y (io- 4 ^ 1 - 2 

1 + 10" 2 x 



and 



y = 




if x < 10~ 4 
if 10" 4 < x < 1 
if x > 1 



(11) 



(12) 



(13) 



Fig. [I] compares the photoelectric heating rate derived 
in this paper for large graphite grains with that derived 
by Bakes & Tielens (Ofl for ISM graphite 



for electron particle densities 10 5 < n c < 10 5 cm 
temperatures 10 < T < 10 000 K, and FUV photon fluxes e 
If) 



r po = 10~ 24 e x n H , (14) 
where the photoelectric efficiency e is approximated by 
0.0487 0.0365 (10- 4 T) 7 



1 



1 + 210- 



(15) 



The difference is due to the larger photoelectric yield of 
small grains. Photoemission is mostly a "volume process" . 
In large grains the escape probability of photoelectrons is 
reduced, because they are created far below the surface. 
There scattering becomes important and they can easily 
be retrapped in the grain. Fig. || reveals the same effect 
for the large silicate grains. 
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Heating by collisional de-excitation of H2. H2 molecules ex- 
cited by Lyman and Werner band absorption decay to vib- 
rotationally excited levels of the ground electronic state 
and heat the gas via collisional de-excitation. We approx- 
imate the complex process using the formula of Tielens & 
Hollenbach (1985) for a single excited pseudovibrational 
level 

r 2 = /h- n(R 2 ) (l.67KT n n(H) VT c - 1000 / t + 
2.33 HT 11 n(H 2 ) \Pf e - 18100 /( T + 120 °) N 



— 3 — 1 

erg cm s 



(16) 



/h* denotes the fixed fraction of vibrationally excited 
molecular hydrogen and E* = 4.166 10~ 12 erg is the ef- 
fective energy of the pseudolevel. We assume in this paper 
fm = 10- 5 . 



Heating due to C ionisation. In the presence of strong 
UV radiation neutral carbon can be photoionized with 
an average yield of approximately 1 eV per photoelectron 
(Black |i~987|) 



1.6022 10' 12 T c n(C) / shicld (C) X x /shicid(H 2 ) 



erg cm 3 s 1 



(20) 



where Xx is the strength of the UV radiation field atten- 
uated by dust. The carbon self-shielding factor 



/shield (C) = exp (-JV(C) oc) 



(21) 



depends on the C column density iV(C) and on the average 
carbon ionization cross section ac- The a ttenu ation by H2 
can be approximated (Kamp & Bertoldi 200C| ) by 



/shicid(H 2 ) = exp ( -T g ° as 25 (iV(H 2 )/10 2 



(22) 



Heating by photodissociation of H 2 . Photodissociation of 
H 2 occurs via line transitions to excited levels followed 
by spontaneous radiative dissociation into two hydrogen 
atoms. The kinetic energy of each H-atom is typically 
0.4 eV (Stephens & Dalgarno 1973] ), leading to an ap- 
proximate heating rate 



r,- 



6.410 



-13 



H 2 



n(H 2 ) 



erg cm 3 s 1 



(17) 



shield 



s 1 denotes the H 2 photodissociation 



rate, which depends on the strength of the UV radiation 
field and on the amount of H 2 self-shielding (Kamp & 
Bertoldi p000| ). 



Heating by H 2 formation on dust. The formation of an H 2 
molecule on the surface of a dust grain releases the binding 
energy of 4.48 eV. Due to the lack of laboratory data, 
we follow the approach of Black & Dalgarno (1976) and 
assume that this energy goes into translation, vibration 
and rotation in equal parts. Hence about 1.5 eV will go 
into kinetic energy and therefore heat the gas 



r 4 = 2.39 10~ 12 i? form n(H) n tot erg cm" 3 s" 



(18) 



-Rform is the temperature dependent H 2 formation rate 



(Kamp & Bertoldi 2000) 



Gas-grain collisions. Gas-grain collisions will act as a gas 
heating source for dust temperatures larger than the gas 
temperature (Burke & Hollenbach 1983| ), 



r 5 = 4.0 10" 12 7ra 2 n tot n dust a T VT (T dust - T) 

erg cm -3 s -1 , (19) 



where ndust an d c^t denote the particle density of dust 
grains and the thermal accomodation coefficient respec- 
tively. A typical value for sili cate and graphite dust is 
ax = 0.3 (Burke & Hollenbach 1983 ). If the gas tempera- 
ture exceeds that of the dust, this rate becomes negative 
and therefore acts as a cooling rate (A5). 



Drift velocity of the dust grains. Assuming typical physical 
conditions for these circumstellar disks, gas and dust are 
momentum coupled (Oilman ^972) and we assume that all 
the momentum gained by the grains from the radiation 
field is transferred to the gas by collisions. Dust grains 
of micron-size are usually removed by radiation pressure 
from the disk on timescales short compared to the life- 
time of the disk (Artymowicz & Clampin 1997; van der 
Blick et al. 1994| ; Backman & Paresce 1993; Artymowicz 
1988). Since observations nevertheless reveal the presence 
of micron-sized dust grains, there must be a process of con- 
tinous replenishment for these small grains, for example 
destructive collisions between much larger dust particles. 

Gas and dust in the disk rotate with slightly subkep- 
lerian speeds v<f, due to the gas pressure (gas component) 
and due to the drag force between dust and gas (dust 
component). Hence stellar gravity is largely balanced by 
centrifugal forces. The evaluation of the drift velocities in 
both directions leads to a coupled system of two-fluid two- 
dimensional hydrodynamical equations and is beyond the 
scope of this paper. Hence we restrict ourselves here to the 
case of maximum drift velocity, that is assuming keplerian 
rotation, and bracket the "true" solution by thrift = an d 
■y drift = v™?ft. The maximum drift velocity can then be de- 
rived from the balance of radiation pressure on the grains 
and dra g forc e due to the grains motion through the gas 
(Tielens |l983|) 



max 
"drift 



2 ( ((/radf 



0.5 



O.o 



cm s 



with the thermal velocity of the gas 



^gas 

and 

frad 



I bkT 



cm s 



2ncr 2 p tot 



cmV 2 



(23) 



(24) 



(25) 
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Here a g and m g are the geometrical cross section and the 
mass of the dust grains, L* and M* are the luminosity 
and mass of the central star, p t ot is the total mass den- 
sity of the gas, and Q GXt is the flux-weighted mean of the 
radiation pressure efficiency 



Qext — 



J °°Qext(^) F{V) dv 

J °°F(v) dv 



(26) 



The rate of viscous gas heating is approximated 
(Tielens |1983[ ) by 

T 7 = 0.5 ptot "dust c g vl Tilt erg cm" 3 s" 1 . (27) 



Cosmic ray heating. Following Clavel et al. ( 197§D The cos- 
mic ray heating rate for a mixture of H, H2 and He is 



r 8 = (1 + i Hc ) 



C(H) n(H) (1.2810" 

„— 1 



-2.4410" n a;H 2 ) 

erg cm " s * , (28) 

where C(H) = 6 10 -18 s _1 is the primary cosmic ray ion- 
isation rate of hydrogen and xu e and Xn 2 are the abun- 
dances of He and molecular hydrogen relative to atomic 
hydrogen. 

2.1.2. Cooling processes 

In the following we give a short description of the cooling 
processes taken into account for the determination of the 
gas temperature. 

Oi cooling. At densities larger than n^ r = 
8.5 10 5 (T/100)-°- 69 cm- 3 , LTE is a good approxi- 
mation for the oxygen level occupation numbers. However 
the densities in the outer parts of these disks will be 
below that critical value. Therefore we have to calculate 
the statistical equilibrium (SE) of oxygen in detail. More 
importantly thermal emission of 50 K dust peaks at 
60 fj,m, where we also find the [O 1] 63.2 /mi fine-structure 
line. Since typical dust temperatures in our disks range 
from 20 to 150 K, this can lead to a strong pumping of 
the fine-structure levels of neutral oxygen. An additional 
component to the infrared radiation field P„ is the cosmic 
microwave background, which is more important for low 
rotational levels of molecules like CO. 

We use an oxygen model atom consisting of the low- 
est three fine-structure levels and we take into account 
the three fine-structure lines at 63.2, 145.6, and 44.0 /itm, 
and collisions with H2, H, and electrons (Jaquet et al. 
1991; Launay & Roueff|l977j Bell et al.p98|). We include 



spontaneous emission as well as absorption and stimulated 
emission due to the IR radiation field. The importance of 
the latter two processes for the statistical equilibri um ca l- 
culations is stressed by Kamp & van Zadelhoff fl200CD . 
Table gives a short overview of the atomic data. The 
last column gives fits to the critical densities above which 
collisional de-excitation becomes important 



A, 



cm 



(29) 



where K ui = e H C* + en 2 C^ + e e -C^ denotes the sum 
of all collisional transitions from level u to all other lev- 
els i. The critical densities tabulated refer to a neutral 
molecular gas with en = 10~ 3 , en 2 = 1.0, and e e - = 10~ 6 . 
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Fig. 3. [0 1] fine structure cooling rates for (a) T gas = 20 K 
and (b) T gas = 100 K; solid line: SE including an IR radi- 
ation field P v = 0.01 5^(80 K) + B„(2.7 K), long dashed 
line: SE including an IR radiation field P y = -B„(80 K) + 
B„(2.7 K), dashed line: SE with P v = 5^(2.7 K), dotted 
line: LTE including P v = 0.01 5^(80 K) + J B„(2.7 K); thick 
lines indicate a negative cooling rate (net heating of the 
gas) 



The total cooling rate is the sum over all three fine- 
structure transitions 



n u (0) (A u i+B u iP u ) - 
ni(O) B lu P v 



erg cm 3 s 1 



(30) 



where A u i is the Einstein probability for spontaneous 
emission, B u i — c 2 /{2hvf ll ) A u i is the Einstein proba- 
bility for stimulated emission, and B\ u = g u /gi B u i the 
Einstein coefficient for absorption. The level population 
numbers n;(0) and n„(0) refer to the lower and upper 
fine-structure level respectively, and hv u i is the energy ab- 
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Fig. 4. Optical depth logr 63 of the [Oi] 63.2 finx line in a 2 M ffi /3 Pictoris disk model: (a) optical depth logr for a 
photon emitted from the stellar surface, (b) optical depth logr for a photon emitted from the disk-midplane 



sorbed from the total infrared radiation field P v or radi- 
ated away per emission. Ai turns into a heating rate as 
soon as is becomes negative. 

Table 1. Line data for the oxygen model atom: 
lower and upper levels of fine-structure line, statistical 
weights, wavelengths, transition probabilities for sponta- 
neous emission, and fits to the critical densities for the 
temperature range 10 K < T gas < 10 3 K 



guga 



desig A [/im] A u i [s J ] n cr [cm 3 ] 



10 3 5 
2 1 

2 1 5 3 P - 3 P 2 



3 Pi- 3 P 2 

3l 



1 3 3 P - 3 Pi 
Po- 



63.2 
145.5 
44.1 



8.87 10" 5 3.0 10 6 T-°- 35 
1.77 10" 5 3.710 5 T-°- 35 
1.28 10" 10 2.7T -0 ' 35 



Fig. g illustrates as an example the influence of an 
80 K IR radiation field on the cooling rate. Assuming 
black body dust emission, Ai is a net heating rate at 
low temperatures. Using a more "realistic" radiation field 
of 0.01 -B„(80 K), which includes the dilution of radia- 
tion due to the low densities in the disk, the heating 
rate approaches the statistical equilibrium value at par- 
ticle densities of 10 7 cm -3 . For higher gas temperatures 
(?g as > Tdust) the dust infrared emission is no longer able 
to heat the gas and LTE is a good approximation for parti- 
cle densities larger than 10 7 cm" 3 . Spontaneous emission 
dominates the energy balance of O I and hence the re- 
sulting cooling rate is only marginally dependent on the 
background radiation field. 

The optical depth of the strongest [Oi] line at 63.2 /im 
can be calculated assuming that all oxygen is in the ground 
level, and using a typical line width of <5vd = 1 krn s" 1 



no = A 



10 7 



A? 9i 



N (O) . 



(31) 



87r 15 5i; d g ' 

An optical depth t\q = 1 is reached for column densities 
of ATq(O) = 3.310 17 cm" 2 . Assuming that all oxygen is 
in atomic form and is neutral throughout the whole disk 
(logeo = —3.49), this converts into a total column den- 
sity of iVtot = I- !° 21 cm" 2 . Using Eq.(4) from Kamp & 



Bcrtoldi ( 2000 ), disks up to a total mass of 2 M ffl are op- 
tically thin for [O i] 63.2 /im radiation. Fig. |] illustrates 
the optical depth of this line in a 2 M® (3 Pictoris disk 
model for photons emitted from the stellar surface (left) 
and from the disk-midplane (right). 

Cn cooling. We take the two level approximation of 
Hollenbach & McKee (198£) for the ground state of sin- 
gle ionized carbon; the energies and statistical weights are 
E = 0.00 erg, E x = 1.27 10" 14 erg and g = 2, g x = 4 
respectively. The 157.7 (im line, corresponding to a tran- 
sition from the first excited level 2P 3 / 2 to the ground level 
2Pi/2, has a cooling rate of 



A-2 = fx n(G n) Aw hfio erg cm 3 s 1 



(32) 



where the fraction of C n in the upper level is calculated 
under the assumption of LTE 



fx = 



gx e" 



-Ei /kT 



go + gx e 



-Et/kT 



(33) 



The critical density for 



LTE 



3.0 10 2 (T/100)- 



0.07 



cm is always lower than the 



minimum density of 10 4 cm 3 in our disk models. The 
Einstein probability for spontaneous emission A\o is 
2.4 10" 6 s" 1 . 

H2 rotational/vibrational line cooling. We use the H2 cool- 
ing function derived by Le Bourlot, Pineau des Forets & 



Flower (1999), which can be applied for gas temperatures 
between 100 and 10 4 K and gas densities ranging from 1 
to 10 8 cm" 3 . Due to the wide energy spacing of the levels, 
rotational excitation of H2 is negligible below 100 K. Le 
Bourlot, Pineau des Forets & Flower take into account 51 
rovibrational energy levels and their collisional excitation 
by H, He and H 2 using recent quantum mechanical cal- 
culations of the cross sections. Following their results the 
cooling depends only weakly on the ortho-to-para ratio. 
We take n(ortho)/n(para) = 1.6, the equilibrium value at 
100 K, to be representative for the low gas temperatures 
in our models (see their Fig. 9). 
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1 (6300 /im) cooling. Line emission from the metastable 
*D level of neutral oxygen efficiently cools the gas at 
very high temperatures. We take the cooling rate from 
Sternberg & Dalgarno (|1989[) 



A 4 = 1.810- 24 n(O) n c e - 2280 °/ T 



erg cm 3 s 1 



(34) 



where n(0) and n e denote the density of neutral oxygen 
and electrons respectively. 

Lya cooling. The cooling by Lya emission becomes effi- 
cient at very high temperatures (Sternberg & Dalgarno 



A 6 = 7.3 10~ 19 n n(H) e 



-118400/T 



erg cm 3 s 1 



(35) 



n(H) and n c denote the density of neutral hydrogen and 
electrons respectively. 



CH cooling. Rotational cooling by CH is calculated using 
t he op tically thin approximation of Hollenbach & McKee 
1979b 



4(fcT) 2 A 



n-tot#o (1 + (n-cr/n-tot) + l-5(n cr /n to t) - 5 ) 
erg s _1 cm 3 if n tot » n CI 



(36) 



kT(l - n(H 2 )/n t ot)crtotWT 



1 



erg s 1 cm 3 



where n cr , A and E are 6.6 10 9 T 3 ±/z cm 3 , 7.710 "s 



("tot/n-cr) + 1.5(n to t/"cr) - 5 
if n tot < n cr 

-1/2 _q 

' cm , 

and 2.76 10 -15 erg respectively for the CH radical, vt 
is the thermal velocity of colliding hydrogen atoms. The 
cooling rate is then 

A 7 = ntot "CH L mt erg cm -3 s" 1 , (37) 



where n CH denotes the particle density of the CH 
molecule. 




4 6 
log n m [cm" 3 ] 




4 6 
log n m [cm" 3 ] 

Fig. 5. CO rotational line cooling rates for (a) T gas = 20 K 
and (b) T gas = 100 K; solid line: SE including an IR radi- 
ation field P y = 0.01 S^(80 K) + B„(2.7 K), long dashed 
line: SE including an IR radiation field P v = -B„(80 K) + 
5^(2.7 K), dashed line: SE with P v = B v {2.1 K), dot- 
ted line: LTE including P u = 0.01 B„(80 K) + 5^(2.7 K), 
dash-dotted line: analytic approximation of McKee et al 



(1982); thick lines indicate a negative cooling rate that 
means a net heating of the gas 



CO cooling For the CO molecule the occupation numbers 
of the rotational lines in the ground-vibrational state dif- 
fer significantly from LTE for densities below 10 7 cm' 3 . 
Therefore we use a CO model molecule to calculate the 
statistical equilibrium of carbon monoxide. It consists of 
26 levels up to J=25 and has 351 CO-H 2 collisional rate 



coefficients (Schinke et al. 1985). The pumping of CO lev- 
els by IR radiation is included in the same way as for 
O I (see respective paragraph) . The molecular line data is 
summarized in Table |^. The last column of this table gives 
the best fits to the critical densities for the CO rotational 
lines using the same assumptions as made for the critical 
densities of O I. 

CO rotational cooling is calculated in the optically thin 
limit by summing over all 25 rotational line transitions 



2H 



m(CO)(i4« +B ij P u ) - 



n i (CO)B Si P v 



erg cm 3 s 1 



(38) 



with j 



1. ni(CO) denotes the occupation number for 



the level J and Aij is the spontaneous emission probability 
for the rotational transition J— >J-1. Bij and Bji are the 
Einstein coefficients for stimulated emission and absorp- 
tion respectively and hvij is the released or absorbed en- 
ergy. For photons travelling along the disk midplane, some 
of the CO lines will become optically thick. Nevertheless 
this does not influence the cooling rate, as the energy is 
radiated isotropically in all directions and most of it can 
therefore escape the flat disk without any reabsorption. 



McKee et al. (1982) determined numerical fits for the 



CO cooling rate for gas temperature larger than 150 K. 
These cooling rates are derived under the following as- 
sumptions: 1. optically thin lines, 2. infrared pumping by 
dust emission is negligible, 3. only the ground vibrational 
state contains a significant population. Assumptions 1 and 
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log^. [cm] logA, [cm] 

Fig. 6. Infrared radiation field for a 2 M® (3 Pictoris disk at four different grid points: a) r — 50 AU, disk midplane, b) 
r = 50 AU, z = 15 AU, c) r = 307 AU, disk midplane, d) r = 307 AU, z = 92 AU. The black dots denote the radiation 
field due to the dust derived from Eq. ([io|) . The upper solid curve shows the local approximation P v (r, z) = B v (Td(r, z)), 
the lower one the cosmic microwave background B v (2.7 K). The dashed curve shows the local approximation scaled 
to the shortest wavelength point 



3 hold also for our calculations, but contrary to them, we 
included pumping by the infrared radiation field. Our cal- 
culations show that pumping due to thermal dust emis- 
sion of the disk and the cosmic microwave background 
can be important for the CO population numbers (Kamp 
& van Zadelhoff p000| ). Comparing our cooling rates with 
IR pumping by the cosmic microwave background ( Fig. 5 
dashed lines) to the numerical fits of McKee et al. ( 1982 ) 
(Fig. 5 dash-dotted lines), they agree reasonable well 
within a factor 2.5/3 for T gas = 100 K/20 K. This is within 
the error which the authors claim for their fits. 



Besides this good agreement between our calcula- 
tions and those of McKee et al. ( 1982| ), Fig. 5 shows 
the effect of the IR radiation field P v (solid line: with 
P v = 0.01 B v {m K) + B v (2.7 K), dashed line: with P v = 
B v {2.7 K), long-dashed line: P v = B u (&0 K)+B„(2.7K)). 
At low particle densities the cooling rate turns into a heat- 
ing rate for the two cases P v = 0.01 £„(80 K) + B v (2.7 K) 
and P v = B u {2.7 K). Only at low temperatures and with a 
much stronger radiation field P v — B„(80 K) + B V (2.7 K) 
we find a net heating of the gas over the whole range of 
densities considered here. For the optically thin disk mod- 



els presented here, Fig. 5 and ^| reveal that the pumping 
of the CO molecule by thermal dust emission is negligi- 
ble for densities larger than 10 3 cm -3 and that LTE can 
overestimate the cooling rates by a factor 3 to 10 even at 
high densities. 



[Ci] fine structure line cooling. Considering the chemical 
structure of the j3 Pictoris disks (Kamp & Bertoldi 2000 ) 



[C i] fine structure line cooling may become important due 
to the presence of neutral carbon. The calculation of the 
cooling rates makes use of the three level atom approx- 
imation from Hollenbach & McKee ( |l98Sf ) Ci( 3 P , 3 Pi, 
3 P2). The respective cooling rates for the three lines at 
609.2, 229.9, and 369.0 /zm are 



A 9/io/n = Aijhvijrii erg cm 



3 s- 1 



(39) 



The spontaneous transitions rates Ay are 7.9 10~ 8 , 
2.0 10~ 14 , and 2.7 10 -4 . For the level population numbers 
we can assume LTE, because the gas densities in these 
disks are well above 2 10 3 cm" 3 , the crit ical d ensity for 
LTE (see Table 8 of Hollenbach & McKee |l989|). 
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where <5d g is the dust-to-gas mass ratio and /9 gra in is the 
material density of the dust grains. The total infrared ra- 
diation field is then the sum of the two components 



Table 2. Molecular line data for the CO model molecule: 
upper level of the rotational lines J — > J — 1, energy of 
the upper level, wavelengths of the lines, Einstein transi- 
tion probabilities for spontaneous emission, and fits to the 

critical densities which hold for 10 K < T gas < 100 K; the P v {r, z) = P v ^{r, z) + B v {2.7 K) 
statistical weight of a level is g = 2J + 1 



(43) 



.1 


E [cm" 1 ] 


A [mm] 


A<j [s 


'] 


n cr [cm 


3 ] 


1 


3.85 


2.601 


7.1710" 


-8 


5.0 10 3 T" 


-0.66 


2 


11.54 


1.300 


6.8710" 


_7 


1.9 10 4 T" 


-0.45 


3 


23.07 


0.867 


2.48 10" 


-6 


4.6 10 4 T" 


-0.35 


4 


38.45 


0.650 


6.09 10" 


-6 


8.3 lO 4 ! 1 " 


-0.28 


5 


57.67 


0.520 


1.22 10" 


_5 


1.3 io 5 r" 


-0.22 


6 


80.74 


0.434 


2.13 10" 


_5 


2.0 10 5 T" 


-0.18 


7 


107.64 


0.372 


3.40 10" 


_5 


2.9 10 5 T 


-0.16 


8 


138.39 


0.325 


5.11 10" 


_5 


4.1 10 5 T" 


-0.14 


9 


172.98 


0.289 


7.2910" 


_5 


5.5 10 5 T" 


-0.12 


10 


211.40 


0.260 


1.00 10" 


-4 


7.2 10 5 T 


-0.11 


11 


253.67 


0.237 


1.33 10" 


-4 


9.2 io 5 r 


-0.10 


12 


299.77 


0.217 


1.73 10" 




1.1 10 6 T 




13 


349.70 


0.200 


2.1910" 


-4 


1.410 6 T 


-0.08 


14 


403.46 


0.186 


2.72 10" 


-4 


i.7io 6 r 


-0.075 


15 


461.05 


0.174 


3.3310" 


-4 


2.0 10 6 T" 


-0.07 


16 


522.48 


0.163 


4.02 10" 


-4 


2.4 10 6 T" 


-0.065 


17 


587.72 


0.153 


4.8010" 


-4 


2.810 6 T 


-0.06 


18 


656.79 


0.145 


5.6610" 


-4 


3.2 10 6 T 


-0.055 


19 


729.68 


0.137 


6.6010" 


-4 


3.710 6 T 


-0.05 


20 


806.38 


0.130 


7.6410" 


-4 


4.2 10 6 T" 


-0.045 


21 


886.90 


0.124 


8.7710" 


-4 


4.710 6 T 


-0.04 


22 


971.23 


0.119 


9.9810" 


-4 


5.310 6 T 


-0.035 


23 


1059.37 


0.113 


1.1310" 


-3 


5.910 6 T 


-0.03 


24 


1151.32 


0.109 


1.2710" 


-3 


6.5 10 6 T 


-0.025 


25 


1247.06 


0.104 


1.42 10" 


-3 


6.8 10 6 





3. The infrared radiation field 



We stressed already in Sect. 2.1 the importance of the 
infrared radiation field for the determination of the sta- 
tistical equilibrium of oxygen and carbon monoxide. This 
radiation field consists of two components: one is the ther- 
mal emission of the large dust grains and the other is the 
2.7 K cosmic microwave background. 

Since the disks are optically thin, we can calculate the 
mean radiation field from the dust grains at any point 
(r,z) in the disk straight forwardly by integrating over the 
whole disk volume 

/ -^-n d (r,z) 

-z jRi y 

B u (T d (r' , z'))r' d<j)dz' dr' , (40) 



Pu,d(r,z) = 



1 

8tt2 



where Z a , i2; and R a define the outer boundary of the disk 
model. The distance from (r, z) to any other point (r', z') 
in the disk is given by 



y = + r ' 2 ~ 2rr' cos <j) + (z' — z) 2 , (41) 

and the dust particle density can be derived from the gas 
density 

3 



rid = <5d g n g /i ibh 



47ra 3 /9 g 



(42) 



The dust particle density is too low to ensure a total 
47r coverage even in the innermost regions of the disk. For 
an oxygen atom sitting in the disk plane the surrounding 
medium is more or less a dark background with a lot of 
bright spots, the dust grains. Hence the local approxima- 
tion P v {r, z) — B v (T d (r, z)) completely fails. Fig. || reveals 
the difference beween these two approaches for four dif- 
ferent points in the disk. While in the inner regions of the 
disk the IR radiation field has at least the same tempera- 
ture as the local gas, this is clearly not the case in the cool 
outer parts (Fig. |^). Moreover the radiation field cannot 
be represented by some kind of mean temperature scaled 
according to the solid angle of dust grains. 

Black body grains with a size of 3 /xm have at 60 /im 
an absorption efficiency of Q a bs = 0.6. The assumption of 
a black body radiation field is no longer valid when the 
dust temperature changes by more than 5% over an op- 
tical depth of 1, which means a total column density of 
7 10 22 cm~ 2 . For the disk model described here, the as- 
sumption of a black body radiation field is valid at 50 AU 
for total disk masses (between 40 and 500 AU) larger 
than 7.4 10 4 M ffi and at 100 AU for masses larger than 
4.2 10 5 Mg. These values are well above the disk masses 
considered here. 

The infrared radiation field depends only on the den- 
sity structure of the disk and on the dust temperature. 
Since the latter is fixed by the assumption of radiative 
equilibrium, it can be calculated in advance and does not 
change during the numerical calculation of the gas chem- 
istry and gas temperature. 



4. The gas temperature: numerical models 

The input parameters listed in Table remain the same 
during the calculations presented here; we assume that 
the grains are composed of silicates. The only parame- 
ters changed are the disk mass and stellar radiation field. 
The models are calculated for two different disk masses, 
namely 2 M ffi and 0.2 M ffi and for the radiation fields of 
two prototype stars (3 Pictoris and Vega, which represent 
the cool and the hot end in the range of A stars. 

The radiation field of (3 Pictoris is strong enough to 
blow out 3 fim silicate grains, since these disks rotate with 
slightly subkeplerian velocity around the star. The mo- 
tion of the dust particles will be much more complicated 
than described in Eq.(p3|) and the velocities will probably 
be somewhat smaller. But in order to bracket reality, we 
have carried out the calculations for two extreme cases: 
1. "maximum drift velocity", 2. Udrift = 0. We will first 
discuss the numerical results including thrift = ^drift anc ^ 
afterwards the results without drift velocity heating. 
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Table 3. Standard model parameters 
Parameter Name Standard value 



gram size a 

grain density Pgrain 

inner radius Ri 

outer radius R 

dust-to-gas mass ratio 5d g 

disk scale height H 

fraction of H2 /h* 
UV dust ext. cross section <ttjv 

carbon abundance ec 

oxygen abundance eo 



3 fim 
3 g cm" 
40 AU 
500 AU 
0.01 
0.15 
10" 5 
2.34 1(T 
1.4 10~ 4 
3.2 10~ 4 



(silicates) 



(H)" 



4.1. (3 Pictoris: ^drift 



..max 
y drift 



Drift velocity heating dominates throughout the whole 
disk independent of the disk mass in the j3 Pictoris mod- 
els. The most important cooling processes for the 2 M e 
model are illustrated in Fig. [?]. Since the stellar radiation 
field is too weak to efficiently dissociate carbon monoxide 
(Kamp & Bertoldi 2000| ), the CO rotational lines are the 
most efficient coolant in the inner parts of the disk. Using 
the laboratory data for freezing out of CO on water ice 
surfaces (Sandford & Allamandola 1990), at Tdust < 50 K 
all CO is incorporated in ice mantles around cold dust 



grains (Kamp & Bertoldi 2000); hence [Oi] fine structure 
cooling dominates in those parts. In the uppermost layers 
of the disk from inside to outside [O 1] , [C 11] and then H 2 
line cooling are the most important cooling processes. By 
lowering the disk mass by a factor ten (see Fig. |8|) the dom- 
inant processes become very similar to those found in the 
upper layers of the more massive disk model. The global 
picture does not change much, except that now [C 11] and 
H2 line cooling become more important compared to [O 1] 
cooling. 



< 




1 



300 
r[AU] 



500 



Fig. 7. The most important cooling processes in the 
(3 Pictoris 2 Me disk model including all described heating 
and cooling processes (the bar displays only the relevant 
processes) 

The resulting temperature structure of both disk mod- 
els is shown in Figs. || and Due to the strong drift ve- 
locity heating, the gas temperature of the disk is always 
larger than the dust temperature, except for a small core, 




Am 



300 
r[AUl 



500 



Fig. 8. The most important cooling processes in the 
j3 Pictoris 0.2 M ffi disk model including all described heat- 
ing and cooling processes (the bar displays only the rele- 
vant processes) 



which is below 10 K at 100 AU in the 2 M© model. In 
both models we find a small hot surface layer in the inner 
parts of the disk. We will show later in Sect. ^ that the 
total mass contained in this hot layer is less than 1% of 
the total disk mass. 




300 
r[AU] 

Fig. 9. Temperature structure of the /3 Pictoris 2 M® disk 
model including all described heating and cooling pro- 
cesses: the grey colors show the gas temperature as noted 
on the scale at the right hand side, while the overlayed 
white contour lines show the dust temperature in steps of 
10 K 
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Fig. 10. Same as §, but for a (3 Pictoris 0.2 M© disk 
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4.2. [3 Pictoris: v dri{t = 

If we neglect the drift velocity heating in the (3 Pictoris 
disk, cosmic ray heating becomes important in large parts 
of the disk (see Fig. |ll|). In the inner regions of the disk 
plane, the gas is heated by pumping of the O I levels. 
Looking at the upper layers, we find that in the inner 
regions heating by H 2 formation dominates, while out- 
wards H 2 dissociation is more important. On the other 
hand Fig shows, that [Cil] fine-structure cooling domi- 
nates in the upper disk layers, where carbon is not bound 
in CO. While in the innermost parts CO rotational line 
cooling contributes largely to the total cooling rate, more 
outwards, where CO is frozen onto the cool dust grains, 
[Oi] fine structure cooling determines the gas tempera- 
ture. 



250 




r[AU] 

Fig. 13. The most important heating processes in the 
(3 Pictoris 0.2 M ffl disk model including all described heat- 
ing and cooling processes except the heating due to the 
drift velocity of grains through the gas (the bar displays 
only the relevant processes) 




r[AU] 

Fig. 11. The most important heating processes in the 
(3 Pictoris 2 M ffi disk model including all described heat- 
ing and cooling processes except the heating due to the 
drift velocity of grains through the gas (the bar displays 
only the relevant processes) 




r[AU] 

Fig. 14. The most important cooling processes in the 
(3 Pictoris 0.2 M e disk model including all described heat- 
ing and cooling processes except the heating due to the 
drift velocity of grains through the gas (the bar displays 
only the relevant processes) 
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Fig. 12. The most important cooling processes in the 
(3 Pictoris 2 M e disk model including all described heat- 
ing and cooling processes except the heating due to the 
drift velocity of grains through the gas (the bar displays 
only the relevant processes) 



The heating/cooling processes for the 0.2 M e model 
are shown in Figs. |l3| and [h| Cosmic ray heating is now 
also important in the uppermost low density disk layers. 
Between the two regions of cosmic ray heating, the gas is 
heated by H2 dissociation, while in the inner region of the 



disk H2 formation is more important. The cooling rates 
show a similar behaviour as for the 2 M 9 model except 
that CO rotational line cooling and [0 1] fine structure 
cooling are now more concentrated towards the disk mid- 
plane. 
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1 H a 40 




- 20 




J 





300 
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Fig. 15. Temperature structure of the (3 Pictoris 2 M ffi 
disk model including all described heating and cooling pro- 
cesses except the heating due to the drift velocity of grains 
through the gas: the grey colors show the gas temperature 
as noted on the scale at the right hand side, while the 
overlayed white contour lines show the dust temperature 
in steps of 10 K 
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100 200 300 400 500 

r[AU] 

Fig. 16. Same as ||, but for a /3 Pictoris 0.2 M© disk 

Without any drift velocity heating, the gas in the disks 
becomes extremely cool, T gas < 60 K (Figs. [l5| and pf ). 
Note that one of the main heating sources, namely heating 
by H2 formation, is not as efficient as in molecular clouds, 
because the dust grains used in our disks are much larger 
than in the ISM, and hence the dust formation rate is 



lower (Kamp & Bertoldi 2000) 



4.3. Vega: u d rift 



7 drift 



Contrary to the (3 Pictoris model, CO is mainly photodis- 
sociated in the Vega disks and hence CO rotational line 
cooling is not important. Instead the main coolants are 
1 and C II. 

Similar to the (3 Pictoris models, drift velocity heat- 
ing is the main energy input for the gas. Fig. 17 shows, 



that [Oi] (dense inner parts) and [C11] 157.7 /xm (outer 
parts) fine structure lines are the most important cooling 
processes of the 2 M© Vega model. 




n 100 - 
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Fig. 17. The most important cooling processes in the Vega 
2 M© disk model including all described heating and cool- 
ing processes (the bar displays only the relevant processes) 

If we reduce the density in the disks by a factor 10 
(Fig. |l8|) , the most important cooling processes remain 
the same. 

Fig. |l9] reveals that the gas temperatures throughout 
most of the 2 M© disk model stay well below 100 K. 
Nevertheless the heating due to the drift velocity of dust 
grains through the gas leads to a hot surface layer, es- 
pecially in the inner parts of the disk. Due to the lower 
particle densities in the 0.2 M© disk model this hot surface 



300 
r[AU] 



500 



Fig. 18. The most important cooling processes in the Vega 
0.2 M© disk model including all described heating and 
cooling processes (the bar displays only the relevant pro- 
cesses) 



layer reaches deeper into the disk, leading to temperatures 
above 100 K even in the disk midplane. Most of the outer 
regions of the disk have nevertheless gas temperatures be- 
low 100 K (Fig. [2(] note the different temperature scale). 
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Fig. 19. Temperature structure of the Vega 2 M© disk 
model including all described heating and cooling pro- 
cesses: the grey colors show the gas temperature as noted 
on the scale at the right hand side, while the overlayed 
white contour lines show the dust temperature in steps of 
10 K 
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Fig. 20. Same as Fig. [Lj, but for a 0.2 M© disk 
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4.4. Vega: u drift = 

If the drift velocity is turned off, H2 dissociation and pho- 
toelectric heating take over (Figs. pl| and |23]) . Throughout 
most of the 2 M ffi disk model H 2 photodissociation is the 
main heating source. Only in the uppermost layers photo- 
electric heating by the micron-sized grains become impor- 
tant. As already pointed out in Sect. 2.1, due to the large 



size of the dust grains, the photoelectric heating rate is 2 
to 3 orders of magnitude lower than in molecular clouds 
with typical ISM grains. The main energy loss for this disk 
model occurs via the [Cn] 157.7 fim fine structure line. 
Only in the densest parts of the disk midplane the [Ol] 
fine structure lines are more efficient in radiating away the 
energy (Fig. |2|). 
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Fig. 21. The most important heating processes in the 
Vega 2 Me disk model including all described heating and 
cooling processes except the heating due to the drift ve- 
locity of grains through the gas (the bar displays only the 
relevant processes) 
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Fig. 22. The most important cooling processes in the Vega 
2 Mq disk model including all described heating and cool- 
ing processes except the heating due to the drift velocity of 
grains through the gas (the bar displays only the relevant 
processes) 



If we now reduce the disk mass by a factor 10, the 
[Cn] 157.7 /im fine structure line is the sole remaining 
cooling process (see Fig. |24|). On the other hand, Fig. ^3] 
shows that the heating takes place via ejection of pho- 
toelectrons from the large dust grains (inner part of the 



modelled disk) and H2 photodissociation (outer parts of 
the disk). 
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Fig. 23. The most important heating processes in the 
Vega 0.2 M ffi disk model including all described heating 
and cooling processes except the heating due to the drift 
velocity of grains through the gas (the bar displays only 
the relevant processes) 
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Fig. 24. The most important cooling processes in the Vega 
0.2 Mq disk model including all described heating and 
cooling processes except the heating due to the drift ve- 
locity of grains through the gas (the bar displays only the 
relevant processes) 

The resulting temperature structure of both disk mod- 
els, the 2 M » a nd the 0.2 M ffi disk model, is shown in 
Figs. and P6[ The hot surface layers, which we found 
for the models including the grain drift velocity, have now 
vanished. The whole disk models have moderate temper- 
atures between 10 and 100 K. Even though the gas tem- 
peratures are in the same range as the dust temperatures, 
they do not show the same structuring in the disk models. 
While the dust temperature depends only on the distance 
from the star, the gas temperature depends on the particle 
densities and chemical structure of the disk. 

5. The gas temperature: analytical approximations 

In the inner core of the 2 M ffi (3 Pictoris disk all carbon is 
in the form of CO, hence radiation loss by the rotational 
lines of CO is the main cooling process. On the other hand 
the main heating process is gas heating due to the drift of 
dust grains through the gas. 
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Fig. 25. Temperature structure of the Vega 2 M© disk 
model including all described heating and cooling pro- 
cesses except the heating due to the drift velocity of grains 
through the gas: the grey colors show the gas temperature 
as noted on the scale at the right hand side, while the 
overlayed white contour lines show the dust temperature 
in steps of 10 K 



250 




100 

: 80 








1 H°" 40 




- 20 




J 





100 



200 300 
r[AU] 



400 



500 



Fig. 26. Same as Fig. ||, but for the Vega 0.2 M© disk 
model 



In Sect. 2.1 we find that the analytical formula of 
McKee et al. ( 1982] ) is a good approximation for the de- 
tailed statistical equilibrium CO cooling rate, even at low 
temperatures. The density is always larger than the criti- 
cal density n cr of McKee et al. ( [1982 ) 

cm~ 3 , (44) 



4rri0.75 



= 1.86 10 4 T 



which is a fit parameter in their formula and therefore 
different from the critical densities derived in this paper 
for each CO rotational line. We equate their Eq.(5.5a) 
modified by a factor of 0.5 (see Fig. 5) to our Eq. (|27|) . 
This gives 



'256.6 



3.710- 

eco 



<5dgPto 



^drift 
^Pgrain 



K 



(45) 



where eco = n-co/n-tot is the CO abundance, Sd g the dust- 
to-gas mass ratio, ptot the total gas density, and a and 
Pgrain are the grain radius and the material density of the 
dust grains respectively. Fig. |27] shows the relative error 
(2gas — Tgas)/2g as of this approximation throughout the 
inner disk regions. 

CO is photodissociated in the Vega disks and further- 
more carbon is mainly ionized. In the presence of drift 
velocity heating, we can thus make the simplifying approx- 
imation A2 = I\, [Cll] fine structure line cooling equals 



Fig. 27. Relative error of the analytical approximation for 
the gas temperature given in Eq. ( fl5| ) in the inner parts of 
the (3 Pictoris 2 M© model: e = (T gas - T gas )/T gas in % 



drift velocity heating for the upper disk layers. For low 
gas temperatures the population number of the upper C 11 
level is approximately 2e _9198 / T . This leads to the ex- 
pression 



Tgas — 



91.98 



70.155 + ln(ap grain ecn) - l n ("totVd rift ) 



K (46) 



for the approximated gas temperature. Fig. |2^ illustrates 
that this approximation can be used throughout huge 
parts of the Vega 2.0 M© disk model. Moreover it holds 
in the cool outer regions of the Vega 0.2 M© disk model. 
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Fig. 28. Relative error of the analytical approximation for 
the gas temperature given in Eq.(Br1) in the Vega 2.0 M© 



model: 



(Tgas — 7gas)/Tgas in % 



If we neglect drift velocity heating in these models, we 
can approximate the gas temperature in the inner parts of 
the 0.2 M© model by simply assuming A2 = T\, [Cll] fine 
structure line cooling equals photoelectric heating. For low 
gas temperatures we derive 



91.98 



18.965 + In ecu - In x 



K 



(47) 



This approximation allows to calculate the disk temper- 
ature throughout the whole Vega 0.2 M© model with an 
accuracy better than about 25% (Fig. G9) . 
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Fig. 29. Relative error of the analytical approximation for 
the gas temperature given i 
model: e = (T gas - T gas )/T g; 



the gas temperature given in Eq.(47) in the Vega 0.2 

in % 



i- gas 



6. The gas chemistry 

In order to investigate the temperature dependence of our 
chemical network, we calculate the same models as dis- 
cussed in Sect. ^ but with the assumption T gas = Ti U st- 

To illustrate the effect of the gas temperature on the 
disk chemistry, we show in Fig. |^ and Fig. |l| the CO 
abundances in the (3 Pictoris 2 M® model with the ap- 
proximation T g as = Tdust and with T gas derived from the 
heating/cooling balance. The difference in CO abundance 
is very small and the CO mass in these models is 1.61 10 -3 
and 1.71 1CT 3 M e respectively. 
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Fig. 31. lognco/^tot in the 2 M e model of (3 Pictoris 
deriving T gas from the heating/cooling balance including 
all relevant processes (udrift 



max^, 
U driftJ 



by H2 lines assuming the disks to consist of purely molec- 
ular hydrogen. This is at least for the (3 Pictoris models 
with the heating due to the drift velocity of dust grains a 
significant percentage of the total disk mass. 

Table 4. Overview of disk masses in a certain temperature 
range for all models denotes heating due to the drift 
velocity of dust grains through the gas); all masses are 
given in M® 
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Fig. 30. lognco/n-tot in the 2 M g 
assuming T gas = Tdust 



model of j3 Pictoris 7. Discussion 



This shows that gas chemistry does not change very 
much for gas temperatures below 300 K. Above this 
threshold reactions with an activation barrier start to be- 
come important and they can change the overall chem- 
istry. Table ^ gives the mass of material in certain tem- 
perature ranges. The mass of material hotter than 300 K 
is negligible in our disk models except in the 0.2 M® Vega 
model with drift velocity heating, where it is 5% of the 
total mass. On the other hand a large amount of the total 
mass can have gas temperatures below 30 K. Note that a 
low gas temperature is a necessary, but not sufficient cri- 
terium for freezing out of molecules. For this process the 
dust temperature plays a crucial role. The mass of ma- 
terial hotter than 100 K is the mass that can be traced 



Throughout the whole paper we discussed the two extreme 
cases for the dust grain drift velocity. But is there not a 
better and more accurate way to derive the drift velocity? 

Klahr & Lin ( 200l| ) solved the equation of motion for 
dust grains rotating in a gaseous disk around a central 
star. For the implementation of the hydrodynamical drag, 
it is assumed that the dust grains are in the Epstein 
regime, which means that their friction time scale is much 
smaller than their orbital timescale. Unfortunately this 
is not fullfilled for the 3 /jm grains in our disk models. 
Nevertheless using their Eq.(17) for the radial drift veloc- 
ity, we end up with drift velocities which are a factor 2-5 
lower than the ones we applied in our models. Since the 
drift velocity enters cubic in the respective heating term, 
this would lead to heating rates that are up to a factor 
125 smaller than the ones assumed in our thrift = ^drfft 
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models. But if we use 2 /im grains instead of 3 /im grains, 
we end up with drift velocities that are even somewhat 
larger than ours. 

To summarize, it is obvious that the heating term due 
to the drift velocity of dust grains through the gas de- 
pends strongly on the exact drift velocity. Unfortunately 
the radial drift velocity depends also strongly on particle 
size, which we cannot fix exactly. Moreover the implemen- 
tation of a more realistic hydrodynamical drag term may 
change the resulting radial drift velocity considerably. In 
any case our two approaches Udrift = «drift and thrift = 
will bracket the real solution. 



8. Conclusion 



The models presented in this paper are a link between 
the PDR models for dark molecular cl ouds (Tielens & 
Hollcnbach 1985; Sternberg & Dalgarno |1989D and mas- 
sive flaring disk models us ed for T Tauri stars and Herbig 
Ae/Be stars (Ai kawa et al. |l996[ Chiang & Goldreich |l997| ; 
D'Alessio et al. 1998 ). The optically thin non-flaring disk 
model s for young A stars presented by Kamp & Bertoldi 
(2000) are extended to self-consistently include the cal- 
culation of the gas temperature from a detailed heat- 
ing/cooling balance. 

The model calculations reveal that gas and dust 
temperatures in these disks are completely different. 
Nevertheless there is still a coupling between them due 
to the IR pumping of O I fine structure and CO rotational 
levels. Since the disks are optically thin, this coupling does 
not result in T gas = Tdust, but instead in level population 
numbers that adjust to the local IR radiation field, which 
is several orders of magnitude below the local Planck func- 
tion. 

Since atomic and molecular transition probabilities are 
known with sufficient accuracy, the main uncertainty in 
the calculation of the gas temperature are the physical 
parameters of the dust grains, like the size and the com- 
position. Except for the bright star (3 Pictoris, where the 
silicate featur e at 10 ^im is clear ly detected (Telesco & 
Knacke 



1991; Aitken et al. 1993), we lack observations 



that could constrain the composition of the dust grains in 
these disks. 

The chemical structure of these disks does not change 
significantly, when we drop the approximation T gas = 
Tdust and calculate the gas temperature from a detailed 
heating/cooling balance. Hence the main results of Kamp 
& Bertoldi ( |2000D concerning the CO abundances in these 
disks do not change. 

The gas temperature, the infrared radiation field and 
the statistical equilibrium calculations are crucial for the 
calculation of emission lines from these models. The re- 
sults of the SE calculations show that level population 
numbers can easily be orders of magnitude different from 
thermodynamical equilibrium. 

We will present emission line profiles and possible trac- 
ers of the gas in these disks in a forthcoming paper. 
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